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Abstract 



We have generalized the scattering-matrix theory of multiple Andreev 
reflections in mesoscopic Josephson junctions to the multi-mode case, and ap- 
plied it to short superconductor/normal metal/superconductor junctions with 
diffusive electron transport. Both the dc and ac current- voltage characteris- 
tics are analyzed for a wide range of bias voltages V. For voltages smaller 
than the supeconducting gap the dc differential conductance of the junction 
diverges as l/y/V. 
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Considerable interest, both in experiments [|T]-[(J and in theory is attracted currently 
to finite-voltage transport properties of mesoscopic Josephson junctions with high electron 
transparency. The mechanism of conduction in such junctions is the process of multiple 
Andreev reflections (MAR) ||T0| , |TT| that takes place at the interfaces of the junction scattering 
region, and bulk superconducting electrodes. Phenomelogically, the MAR processes manifest 
themselves in the so-called "subharmonic gap structure" (SGS) in the current-voltage (IV) 
characteristics of the junction: current singularities at voltages V n = 2A/en, n = 1,2,..., 
where A is the superconducting energy gap of the junction electrodes. Recently developed 
quantitative description of MAR in junctions with arbitrary electron transparency D H 



is based on the usual scattering approach for Bogolyubov-de Gennes equations f!2Hl4 
combined with the idea of acceleration of quasiparticles in the junction by non- vanishing 
bias voltage. Point contacts fabricated with the controllable break junction technique fl,|6|.|T5 



allow for accurate comparison between theory and experiment. In particular, the most recent 
experiment || found that the I — V characteristics of aluminum point contacts can be well 
explained by theory for the whole range of the contact transparencies D, from the tunnel 
junction limit D — > all the way to the ballistic contacts with D — > 1. 

So far the scattering theory of MAR has been formulated only for junctions with one 
propagating electron mode. The aim of this work is to extend the theory to the multi-mode 
case and apply it to diffusive SNS junctions. We show that if the scattering matrix S of the 
junction does not depend on energy on the scale of A, the time-dependent current in the 
junction can be represented as a sum of independent contributions from individual trans- 
verse modes. Therefore, the current depends only on the distribution of the transmission 
probabilities D m of these modes. 

We assume that the junction length L (Fig. 1) is smaller than the coherence length £ of the 
superconducting electrodes of the junction and inelastic scattering length l in : L <C £, l in . 
In this case the transport properties of the junction are determined by the interplay of 
scattering inside the junction characterized by the scattering matrix S, and the Andreev 
reflection with amplitude a(e) at the two interfaces between the junction scattering region 



and bulk superconducting electrodes which are in equilibrium, 



1 £-sign( £ )(£ 2 -A 2 ) 1/2 , |e|>A, 
a(e) = -r\ (1) 

I e - z(A 2 - e 2 ) 1 / 2 , | e |< A. 



Here e is the quasiparticle energy relative to the Fermi level of the electrode. The scattering 
matrix S is a unitary and symmetric matrix 2N x 2N, where N is the number of propagating 
transverse modes supported by the junction, and can be written in terms of reflection and 
transmission N x N matrices r, t: 



S 



I r t \ 



(2) 



where t' = t T , r' = — (t*)~ 1 r t t, and ttf + rr" 1 " = I. 

Due to acceleration of quasiparticles by the applied bias voltage, an electron with energy 
e emerging from one of the electrodes generates electron and hole states at energies e + 2neV 
with arbitrary n. Thus, the electron and hole wavefunctions in regions I and II of the junction 
(Fig. 1) can be written as follows: 

i>A = E n [(a2nA n + J5 n0 )e ikx + Bn e- ikx }e-^ +2neV W h , 
A = E n [A n e ikx + a 2n B n e- ikx }e-^ +2ney W h , 

, , 4>ei = En[C n e ikx + ain^Dne-^e-^ 2 ^^ , 
(II) (4) 
A = E n [a2n+iC n e ikx + J D n e-^] e - i ( £ +( 2 «+ 1 )^)*/ ?i , 

where a n = a(e + neV). In these equations we took into account that the amplitudes 
of electron and hole wavefunctions are related via Andreev reflection. Furthermore, we 
neglected variations of the quasiparticle momentum k with energy assuming that the Fermi 
energy in the electrodes is much larger than A. The quasiparticle energies in regions I and 
II are measured relative to the Fermi level in the left and right electrode, respectively. 

The amplitudes of electron and hole wavefunctions have transverse mode index m not 
shown in eqs. (|]) and (f|), e.g., A n = {A njm }, in = 1, N. The source term J describes an 



electron generated in the jth transverse mode by a quasiparticle incident on the constriction 
from the left superconductor: J(e) = (1— \ do\ 2 ) 1 ^ 2 5 m j. The current I(t) in the constriction 
can be calculated in terms of the wavefunction amplitudes. Equation ([3]) and (f|) imply 
that the current oscillates with the Josephson frequency ujj = 2eV/h and can be expanded 
in Fourier components: I(t) = J2 k Ik£ lkuJj ■ Summing the contributions to the current from 
electrons incident both from the left and right superconductors at different energies e we 
obtain the Fourier components Ik'- 

h = — % r c/etanh{ 7 ^}Tr[(JJ t 4o + a* 2k JA\ + a__ 2fc A_ fc J t + 



+ £(! + 



fl— >oo 



(5) 



where Tr is taken over the transverse modes. 

Amplitudes of the wavefunctions (f|) and (|4|) are related via the scattering matrix S. 
Taking into account that the scattering matrix for the holes is the time-reversal of the 
electron scattering matrix 5", we can write: 
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Eliminating A n between eq. (|6]) and inverse of eq. (^), we find a relation between the ampli- 
tudes B n and D n . Combining this relation with the expression for D n in terms of B n that 
follows from the inverse of eq. (^j and eq. (0) we arrive at the following recurrence relation 
for B„: 



B n 



+i 



h 2n+l 



'2)i 



2n+l 
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Since the hermitian matrix ttf can always be diagonalized by an appropriate unitary trans- 
formation U, the recurrence relation (§) implies that the structure of the amplitudes B n as 
vectors in the transverse-mode space is: 

B n = U^f n (D)UrJ, (9) 

where D = Utt'W is the diagonal matrix of transmission probabilities D m , m = 1,...,N. 
The functions f n {D) are determined by the solution of the recurrence relation (^) with the 
diagonalized transmission matrix tv . 

Equation (||) shows that the contribution of the amplitudes B n to the currents (|5]) can 
be written as 

Tr[B n Bl] = (1- | a | 2 )Tr[/„(L>)/*,(L>)(l - D)) , 

i.e., it can be represented as a sum of independent contributions from different transverse 
modes with the transparencies D m . Following similar steps we can show that the same is 
true for the amplitudes A n . Therefore, the Fourier components (|5]) of the total current can 
be written as sums of independent contributions from individual transverse modes: 

I h = Y,h{p m ), (10) 

m 

where the contribution of one (spin-degenerate) mode is: 



h{D) = \ 



eVD5 k0 - J detanh{^}(l - \a G \ 2 )(a* 2k Al + a_ 2 kA- h 



+ + a 2nO*2(n+k)){AnA* n+k - B n B* n+k )) 



(11) 



with the integral over e taken in large symmetric limits. The amplitudes B n in this equation 
are determined by the recurrence relation which follows directly from eq. (§): 



2 2 

Ogn+gOgrt+l p r n , a 2n+l , a 2n \ , 1 „2]r> 

L>- g ~ 2 1 2 > + 1 ~ a 2n\ B n + 

1 ~~ a 2n+l 1 — a 2n+l 1 — a 2n-l 



+ D a 2n a 2 n-i B _ R i/2 S R ^i- D . (12) 
1 — ai 



b 2n-l 



Instead of using similar independent recurrence relation for A n , it is more convenient to 
determine these coefficients from an equivalent relation that can be obtained from a single- 
mode version of eqs. (|6|) and ([7]) ||: 

- B n a 2n +i) + ai<5„ , (13) 

Equations (0) - (|i~3"D completely determine the time-dependent current in a short con- 
striction with arbitrary distribution of transmission probabilities. In particular, we can 
use them to calculate the current in a short disordered SNS junction with large number 
of transverse modes N ^> 1 and diffusive electron transport in the N region. In this case, 
the distribution of transmission probabilities is quasicontinuous, and is characterized by the 
density function p(D) (see, e.g., |L7j and references therein): 

where G is the normal-state conductance of the N region. 

Figure 2 shows the results of the numerical calculations of the dc current- volt age (IV) 
characteristics and differential conductance of the short SNS Josephson junction based on 
eqs. fllTD - (|TJ|). We see that the IV characteristics have all qualitative features of highly 
transparent Josephson junctions: subgap current singularities at eV = 2A/n and excess 
current I ex at eV 3> 2A. It is instructive to compare quantitatively these features to 
those in the IV characteristics of a single- mode Josephson junctions plotted in ||. Such 
a comparison shows that the magnitude of the excess current in the SNS junction, as well 
as the overall level of current in the sub-gap region correspond approximately to a single- 
mode junction with large transparency D ~ 0.8. At the same time, the subharmonic gap 
structure and the gap feature at eV = 2A are much more pronounced than in a single- 
mode junction of this transparency. The amplitude of the oscillations of the differential 
conductance corresponds roughly to the junction with D ~ 0.4 (although this comparison is 
not very accurate because of the different shapes of the curves). This "discrepancy" reflects 



the two-peak structure of the transparency distribution ([14]) of the diffusive conductor: the 



abundance of nearly ballistic modes leads to large excess and subgap currents, while the 
peak at low transparencies determines the SGS features. 

At large and small bias voltages the time-dependent current through the constriction 
I ((f), where if = 2eVt/h, can be found analytically. At large voltages, eV ^> A, the 
probability of MAR cycles decreases rapidly with the number of Andreev reflections in 
them. This implies that the higher-order harmonics Ik of the current decrease rapidly with 
increasing k, and we can limit ourselves to the first harmonic: 

I ((p) = I(V) + 2Re/i cos(tp) + 2ImJi sin(^) 

Truncating then the recurrence relations (|T2| ) and (|1^) at the coefficients B±i and Ao,i we 
can solve them explicitly and find the current from eq. (pTTj) . For a single mode at T <C A 
we get: 

'00 - - % 



eV H — — (1 =- ln( -=)) 



R 2VR(1 + R) 1-VR 2eV 

Im/ '00-aFTTfl- 05) 



Expression for the excess current (second term in I(V)) was obtained before in [|T~8| , [T9 
Asymptotics of the ac components of the current agree both with the known results for the 
tunnel junction limit (D — > 0) and ballistic junction (D — > 1) ||. 

Averaging eqs. (p~5[ ) with the distribution (0) we obtain the large-voltage asymptote of 
the current in the SNS junction: 

nV) = G[V + ^- !)-£)], (16) 

The second term in equation for /(V) represents the excess current and was first found in 
PI by the quasiclassical Green's function method. It can be checked that the asymptote 
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(|16|) agrees well with the numerically calculated zero-temperature IV characteristic shown 
in Fig. 2a. 

Analytical results at low voltages, eV <C A can be obtained using the understanding || 
that the small voltage V drives the Landau-Zener transitions between the Andreev-bound 
states of the modes with small reflection coefficients R <C 1. Averaging the nonequilibrium 
voltage-induced contribution to the current (eq. (11) of Ref. ||) with the distribution ( |l4l) we 
find the nonequilibrium part of the dynamic current-phase relation of a short SNS junction 
at eV < A: 

«GA,(A)|0. <><„<», 



'eV/A sin(v9/2) , vr < <p < 2tt , 
where 77(A) = tanh (A/2T). (Equilibrium part of the current-phase relation has been found 
before in pUfl .) From this equation we can find the voltage dependence of the Fourier harmon- 
ics of the time-dependent current at low voltages. The amplitudes of the first harmonics 
calculated numerically for arbitrary voltages are shown in Fig. 3. The curves agree with 
the high-voltage fllED and low-voltage fll7D asymptotics and exhibit the SGS singularities 
at intermediate voltages. In general, the curves look qualitatively similar to those for the 
single-mode junction with intermediate transparency D (see Figs. 2b, c in Ref. |§). 

Equation ( |P7D implies that the dc differential conductance of the junction has the square- 
root singularity at V — > 0. 

dl /7(A), 



G{V) = — = ^YG^A/eV. (18) 

This singularity reflects directly the high-transparency peak in the distribution (|14D and is 
not unique to the diffusive SNS junctions. A junction with the strongly disordered tunnel 
barrier should exhibit the same 1/yV singularity with the prefactor 1/2 in eq. fll8"|) 
replaced with 2/tt. Physically, the origin of this conductance singularity is overheating of 
electrons in the junction by MAR. Electrons with energies inside the energy gap traverse the 
junction many times and as a result are accelerated to energies much larger the eV. This 
means that the effective voltage drop across the junction is much larger than V, leading 

8 



to increased conductance. This mechanism of conductance enhancement is qualitatively 
similar to the so-called "stimulation of superconductivity" |22[ (which is one of the plausible 
explanations of the zero-bias conductance singularities |23| in long semiconductor Josephson 
junctions), although quantitatively the phenomena are quite different. The fact that the 
singularity flTH| ) is caused by electron overheating implies that at very low voltages it should 
be regularized by any mechanism of inelastic scattering. Nevertheless, in junctions shorter 
than inelastic scattering length U ni there should be a voltage range where the conductance 
follows the 1/yV behavior. 

In summary, we have developed a theory of multiple Andreev reflections in multi-mode 
Josephson junctions and applied it to the diffusive SNS junctions and disordered tunnel 
barriers. The hallmark of the MAR processes in these systems is the zero-bias 1 / VV singu- 
larity of the dc differential conductance. We have also calculated the low- and high-voltage 
asymptotes of the ac components of the time-dependent current in the SNS junctions. 

We would like to thank K. Likharev for fruitful discussions. This work was supported by 
DOD URI through AFOSR Grant # F49620-95-I-0415 and by ONR grant # N00014-95-1- 
0762. 
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FIGURES 
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FIG. 1. Schematic diagram of the mesoscopic disordered Josephson junction. I and II denote 
the portions of the contact region separated by the scattering region (hatched) where the motion 
of the quasiparticles is diffusive. 
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FIG. 2. The quasiparticle current (a) and the differential conductance (b) versus voltage at 
various temperatures, for a disordered SNS junction. From top to bottom T = 0, 1, 2, 3A. At low 
voltages, the dc current (a) has a square root dependence on voltage, while at high voltages exhibits 
excess current. The conductance (b) posseses subharmonic singularities, diverges at low voltages, 
and at high voltages asympotically tends to the normal state conductance of the disordered region. 
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FIG. 3. The cosine (a) and sine (b) part of the first Fourier component of the ac current. 
From top to bottom (b) and from bottom to top (a) T = 0, 1,2, 3A. The slope at small voltages 
diverges as (A/V) 1 ^ 2 . See text for discussion about the high voltage behavior. 
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